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Abstract. It has been suggested that a continuous low 
star formation rate has been the dominant regime in 
IZw 18 and in dwarf galaxies for the lifetime of these ob- 



reached a particular place. Indeed, in hierarchical clus- 
tering theories these galaxies are the building blocks of 



jects ( pegrand et al. 199E ). Here, we discuss and model 



various star-forming histories for IZw 18. Particularly, we 
show that if the metallicity observed in IZw 18 results 
from starburst events only, the observed colors constrain 
the fraction of the metals ejected from the galaxy to be 
less than 50-70%. We demonstrate that the continuous 
star formation scenario reproduces the observed parame- 
ters of IZw 18. A continuous star formation rate (SFR) 
of about 10 -4 Moyr^ 1 during 14 Gyr reproduces pre- 
cisely the observed abundances. This SFR is comparable 
with the lowest SFR observed in low surface brightness 
galaxies ( Van Zee et al. 1997c ). Generalized to all galax- 
ies, the low continuous SFR scenario accounts for various 
facts: the presence of star formation in quiescent dwarfs 
and LSBG, the metallicity increase with time in the most 
undcrabundant DLA systems, and the metal content ex- 
trapolations to the outskirts of spiral galaxies. Also the 
apparent absence of galaxies with a metallicity lower than 



larger systems by merging (Kauffmann et al. 1993; Pas 
carelle et al. 1996| ; Lowenthal et al. 1997). Moreover, as 



primeval galaxies may undergo rapid and strong star for- 



mation events (Partridge & Peebles 1967), nearby dwarf 



starburst galaxies or Blue Compact Galaxies (BCDG) of 
low metallicity can also be considered as their local coun- 
terparts. Therefore the study of low redshift starbursts is 
of major interest for our understanding of galaxies forma- 
tion and evolution. 

As BCDG presently undergo a strong star formation 
(which cannot be maintained during a long time), but 
generally present a low metallicity indicating a low level 
of evolution, Searle & Sargent (1972) have proposed that 
these systems are young in the sense that they are form- 
ing stars for the first time. An alternative is that they 
have formed stars during strong starburst events separated 
by long quiescent periods. However, most dwarf starburst 
galaxies show an old underlying population indicating that 



they have also formed stars in the past (Thuan 1983: Dou- 



IZw tt , the apparent absence of HI clouds without op- blicr 1998). Thus they are not "young". 



tical counterparts, and the homogeneity of abundances 
in dwarfs galaxies are natural outcomes of the scenario. 
This implies that, even if starbursts are strong and im- 
portant events in the life of galaxies, their more subdued 
but continuous star formation regime cannot be ignored 
when accounting for their chemical evolution. 

Key words: Galaxies - Galaxies: ISM - Galaxies: enrich- 
ment of ISM - Galaxies: - Galaxies: IZwl8 - ISM:Outflows 



1. Introduction 

A challenge of modern astrophysics is the understanding 
of galaxies formation and evolution — In this exploration 



Among starbursts, IZw 18, as the lowest metallicity 
galaxy known locally, could be considered as the best can- 
didate for a truly "young" galaxy. However, recent studies 
have shown that even this object is not forming stars for 
the first time. Color magnitude diagrams have revealed 
the presence of stars o lder than 1 Gyr ( Aloisi et al. 1999 ). 
Legrand et al. (1999 ) have shown that the extreme ho- 
mogeneity of abunda nces throughout the galaxy (see also 
Van Zee et al. 1998| ) cannot be explained by the metals 
ejected from the massive stars formed in the current burst 
(see also Tenorio-Tagle 1996), thus indicating previous 
star formation. Then we need to constrain this previous 
star formation and specify its nature. 

It is generally accepted that the enrichment of the ISM 



arises by burst phases. In the case of IZw 18, Kunth et al 



i j r T- i • u • i (19951) have shown that one single burst with intensity 
low-mass dwarls and irregular galaxies have progressively — i -V . 6 J 
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comparable to the present one is sufficient to reproduce 
the observed abundances. However metals ejected by mas- 
sive stars coul d escape the galaxy, if its t otal mass is lower 
than 10 8 M Q (|Mac Low fc Ferrara 19991 ). The metallicity 
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is then no longer a measure of the number of bursts. On 
the other hand, if the total mass of the galaxy amounts 
M Q the metals are likely to be retained ( |Silich fc" 



10 



Tenorii PTaglc 1998| ). As the total mass of IZw 18 is likely 



to lie between these values ( |Viallefond et al. 1987 ; Van 
Zee et al. 1998]) , the escape of a fraction or of the totality 



of the newly synthesized metals during a burst cannot be 
excluded. In such a case several bursts are needed to ac- 
count for the observed metallicity, their number depending 
on the fraction of metals leaving the galaxy. Nevertheless, 
even if metals escape, stars are likely to remain bound and 
for an increasingly larger number of bursts, the old under- 
lying stellar population will appear progressively redder. 
Thus the number of previous bursts is limited, considering 
the extremely blue color of BCDG. 

Between starburst events, BCDG are likely to appear 
as Low Surface Brig htness Galaxies (LSB G). However, 
studies of the latter ( Van Zee et al. 1997c ) showed that 
despite their low gas density, they do not have a zero star 
formation rate (SFR). LSBG inde ed present a low and 
possibly continuous SFR. This led Legrand et al. (1999 ) 
to propose that a continuous low SFR over a Hubble time 
as responsible for the observed metallicity level in the most 
metal-poor objects like IZw 18. 

Several studies of the past star formation history of 
BCDG, and specifically IZw 18, have been carried out. 



Most have dealt with their chemical evolution (Chiosi fe 



Matte red 1982^ |Carigi et al. 1995| ; [Kunth et al. 1995|) 



or with their spectrophotomctric properties (Mas-Hesse 



& Kurth 1991 



Stasinska & 



Lcitherer fc Heckman 1995 
Leithei er~1996| ; |Cervino fc Mas- Hesse 1994| ; [Mas- Hesse fc 
Kunth 1999), but rarely with both. Moreover, solely the 



influence of bursts has been studied up to now, and the low 
continuous SFR of inter-burst phases has been ignored. 

We used a spectrophotomctric model coupled with 
chemical evolution in order to constrain both the abun- 
dances and the colors of the galaxies. We also investigated 
the effect of a continuous and low star formation regime. 
A preliminary study (Legrand & Kunth 1998) showed that 
this scenario is plausible. Here we present detailed calcu- 
lations, results and their implications. The model and the 
observational data used are described in section ||. The 
different models, including the investigation of mass loss 
effect and the continuous star formation rate model are 
presented in section [| Consequences and generalization 
of the continuous SFR hypothesis are discussed in section 



2. Modelling the star formation history in IZw 18 

In order to investigate the star formation history of 
IZw 18, we used the spectrophotometric model coupled 
with the chemical evolution program "STARDUST" de- 
scribed by Devriendt et al. (1999). The advantage of this 



model is that both the metallicity and the spectral prop- 
erties of a galaxy are monitored through time. 



2.1. The model 

The main features of the model are the following: 

— A normalized 1 Mq of baryonic matter galaxy is con- 
sidered. 

— The SFR and the IMF are fixed and used to evaluate 
at each time the number of stars of all masses formed. 

— The stellar lifetimes are taken into account, i.e., no 
instantaneous recycling approximation is used. Metals 
ejected (C,0, Fe and the total metallicity) are calcu- 
lated at each time step as well as the number of stars of 
each mass. The chemical and spectroscopic evolution 
is followed in time. 

— A fraction of the metals produced by the massive stars 
(M > 9Mq) can be expulsed from the galaxy and do 
not contribute to the enrichment. 

— The fraction of the produced metals remaining into the 
galaxy is assumed to be immediately and uniformly 
mixed with the interstellar medium. We must keep in 
mind that there may be a time delay between their 
production and their visibility. 

— The newly formed stars have the metallicity of the gas 
at the time of their birth. 

— The spectrum as a function of time is computed by 
summing the number of stars multiplied by their in- 
dividual spectra. The nebular emission is not included 
in the model. 

— The model uses the evolutionary tracks from the 
Geneva group ( [Schaller et al. 1992; pharbonncl et al 



1996). The yields are from Macdcr (1992[) for the mas 



sive stars (M > 9M Q ) and from Renzini fc Voli (1981) 



for the lower mass stars. All the metals produced are 
ejected (Case A of |Maeder 19921 ) . 
The stellar output spectra is com puted using the 
stellar libraries from Kurucz (1992 ) supplemented by 
Bessel et al. (1989|, p,99l[ ) for M Giants, and |Brctt 



1995|) for M dwarfs. 



We used a typical IMF described as a power law in the 
mass range 0.1-120 M©. 

4>{m) = a.m~ x (1) 



A constant index x of 1.35 was used (Salpeter 1955) 
We now have some indications that the IMF may flat- 
ten at low masses, maybe below 0.3 M Q ( [Elmegr eer 



1999| ; |5calo 1998j ). As the stars in this range (0.1-0.3 
M©) do not contribute significantly to the enrichment 
of the ISM, nor to the colors, this will only act on the 
normalisation of the SFR in the sense that forming less 
low mass stars will decrease the total SFR requested to 
reproduce the observed abundances. However, as the 
SFR quoted b y |Van Zee et al. (1997b| fl) and repro- 
duced in table 111 are computed using a Salpeter IMF 
down to 0.1 M we used this value in the model in or- 
der to compare our results with these previous studies. 
Finally, their upper mass limit is 100 M Q (against 120 
M in our models). However, the upper mass limit of 
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the IMF do not affect strongly the derivation of the 
SFR but can modify the abundances. Thus using a 
Salpeter IMF ranging from 0.1 M to 120 M© ap- 
pears as a good compromise to study abundances and 
compare SFR with previous studies. 
Two regimes of star formation have been investigated: 

— A continuous star formation during which the SFR 
is low and directly proportional to the total mass 
of available gas. 

— A burst of star formation during which all the stars 
are formed in a rather short time. 



2.2. Comparison of the model with IZw 18 

As the model is normalized to one solar mass of gas, we 
had to multiply the parameters by the mass of IZw 18 in 
order to compare our results with the observations. How- 
ever, this normalization appears only in the value of the 
SFR and in the absolute magnitude predictions; the colors 
reported are independent of the choice of the mass of the 
galaxy. 

The initial mass of gas in IZw 18 must lie between 
6.9 10 7 and 810 s M©, which are respectively the mass 



of HI and the dynamical mass measured by Lequeux & 
Viallefond (1980| ). However, if only the main component 
is considered, the mass must be taken between 2.6 10 7 
(HI) and 2.6 10 s (dynamical mass) as measured by Van 



Zee et al. (1998). The initial mass of gas in IZw 18, in 



the absence of infall, was higher than the mass of HI 
because of the presence of stars and perhaps molecular 
H2 ( Lequeux fc Viallefond 1980| ). As dark matter can 
represent a non negligible fraction of the total dynamical 
mass, the mass of (baryonic) gas can be lower than the 
dynamical mass. We adopted a value of 10 8 M© for the 
initial mass of gas in IZw 18. 

The model produces both the abundances and the 
spectra. We used the spectrum to derive the expected 
colors in (U-B), (B-V) and (V-K). We have adopted for 
comparison with the observed abundance values reported 



by Garnctt et al. (1997) for C and |Skillman fc Kcnnicutt 
(1993| ) for O. Mo st of the pub l ished colors for IZw 18 



are relatively old (Huchra 1977; Thuan 1983), an d have 



not been corrected for the nebular contribution. Balzer 



(1998 ) has recently measured the colors of IZw 18 and 



corrected for the nebular contribution. As our model does 
not include the nebular contribution, we adopted Salzer's 
values for comparison, i.e., (U — B) = —0 .88 ± 0.06 and 
(B-V) = -0.03 ± 0.04; frhuan (1983[) has estimated 



that the flux measured in the IR was mainly of stellar in 
origin. We thus adopted his value for (V-K)=0.57± 0.23. 



Finally, the model used have been compared by pe- 

i.e. 



vriendt et al. (1999) with two similar mo dels, i.e., PE- 
GASE (|Fioc fe Rocca-Volmerange 19971 ) and GISSEL 
( Bruzual A. fc Chariot 1993 ), and no differences larger 



than 0.1 magnitude were found; this was considered the 
intrinsic uncertainty of the modeling process. 



3. Results of the modelisation 

3.1. Enrichment by one previous burst 

It is generally admitted that the starburst events are the 
main contributors to the enrichment of the ISM. We thus 
used the model to evaluate the characteristics of a sin- 
gle burst to reproduce the observed oxygen abundance in 
IZw 18. Taking the uncertainties into account, we found, 



like previous studies (for example Kunth et al. 1995), that 
the present day abundances can be reproduced by a sin- 
gle burst, previous to the current one, with a SFR of 
0.065 Mq yr~ x during 20 Myr. Moreover, the contribution 
of this old underlying population is too faint to modify 
significantly the colors of the galaxy which are currently 
dominated by the newly formed massive stars. This model 
can reproduce all the observations. 

However, some simple arguments can rule out this 
model. Indeed, we assumed that between bursts, the SFR 
is equal to zero, which is certainly wrong. We will demon- 
strate below that even a low but continuous SFR between 
bursts, as observed in LSBG, is likely to produce signif- 
icant enrichment. Moreover, the kinetic energy liberated 
in such a burst is high (about 10 40 erg s" 1 using the mod- 
els of pervino 199"8| ). [Mac Low fc Ferrara (1999| ) have 
suggested that for galaxies with masses comparable with 
that of IZw 18, such an energy is likely to eject out of the 
galaxy all the metals formed by massive stars. If true, it 
means that bursts are unlikely to enrich the ISM by much! 
We thus have investigated the effect of the loss of newly 
synthesized elements by galactic winds. As a less extreme 
hypothesis, we assumed that only a fraction of the SN 
ejecta (and not the totality) leaves the galaxy. 

3.2. The effect of metal loss 

The possibility that the energy released by the SN could 
eject their products out of the galaxy was proposed ear- 
lier by Russell et al. (1988| ). However, intermediate mass 
stars, evolving more slowly, will eject their metals after 
the SN explosion of the most massive stars. Since the ki- 
netic energy released, mainly in stellar winds, is lower than 
for the massive stars, their metal products should be re- 
tained. This will result in a low effective enrichment in 
oxygen (main product of massive stars), but in a rela- 
tively normal enrichment in carbon (mostly produced in 
intermediate mass stars). This hypothesis of "differential 
galactic wind" seems to be necessary to reproduce the 
abundance measurements in some but not all the galaxies 



( [Marconi et al. 1994| ; |Tosi 1998| ). 

If a fraction of metals escapes from the galaxy during 
a burst, the number of bursts necessary to reach the ob- 
served abundance in IZw 18 will be larger. We thus ran a 
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Fig. 1. Time evolution of the gas fraction, oxygen and car- 
bon abundances for recurent bursts with a Salpeter IMF, 
duration of 20 Myr and SFR=0.065 M Q /yr. The dots rep- 
resent the measured abundances. 
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Fig. 2. Colors evolution (U-B, B-V, V-K) for recurent 
bursts with a Salpeter IMF, duration of 20 Myr and 
SFR=0.065 M Q /yr. The plots include the additional 0.1 
uncertainty (from the comparison between the models) 
with the observational errors to define a region (between 
dashed lines) of compatibility between the model and the 
observations. The solid lines represent the prediction of 
the model. If the solid line goes into the region between 
dashed lines, it means that the model does reproduce the 
observations within the error bars. 



model in which 80% of the metals produced by stars more 
massive than 9 M Q left the galaxy and did not contribute 
to its chemical enrichment. Assuming the same parame- 
ters for the bursts as previously, five bursts were necessary 
to reach the oxygen abundance level seen in IZw 18. We 
thus assumed recurrent bursts occurring every 3 Gyr. The 
results of this model are shown in figures g and |. 

Figure [l] shows that if the oxygen abundance is re- 
produced after 5 bursts, the differential winds hypothe- 
sis results in an overproduction of carbon. Moreover, if 
the carbon a bundance in IZw 18 is lower than th e mea- 
surements of Sarnctt et al. (1997 ) , as suggested by Lzotov| 
(1999))] the discrepancy is even larger and completely rules 



out this model. However, the uncertainties on the yields 
remain large (see for example Prantzos 1998, 1999|). For 



example, the detection of WR in IZw 18 (Legrand et al 



1997) can imply that the mass loss rate of massive stars 
are twice the standard ones and the metals produced by 
the intermediate mass stars may also leave the galaxy. We 
thus think that this argument alone is not strong enough 
to invalidate definitively this scenario. 

On the other hand, figure § shows that after four 
bursts, the expected colors (essentially V-K) do not cor- 
respond with those observed. This is due to the fact that 
the old population remaining from the previous star for- 
mation events contributes more and more and reddens the 
colors. This constraint is strong, since it is difficult to ig- 
nore old population. Of course, the constraint from the ob- 
served colors is not really the number of previous bursts 
(because it depends on their strength) but the ratio of 
young stars (formed in the current burst) versus old stars 
(remaining from all the previous star formation events). 
We thus ran another model with recurent bursts of inten- 
sity comparable to the current one ( Mas-Hesse fc Kunth| 
1999) every 1.5 Gyr. This model shows, like in figure [2| 



that the observed colors become incompatible with obser- 
vations after 6 of these bursts. It thus appears that the 
total mass of stars ever formed previously to the current 
burst cannot be larger than 6 times the total mass of stars 
involved in the current burst. Only 2 or 3 of these bursts 
produce enought metals to account for the observed abun- 
dances (if all the metals remain into the galaxy). Thus if 
the present day metallicity results from previous star for- 
mation events with mass loss, this constrains the fraction 
of metals lost by the galaxy to be lower than 50-70%. For 
the same reason, this rules out models inferring a large 
number of bursts in which most of the metals leave the 
galaxy or produce a hot metal-rich halo as suggested by 
|Pantclaki fc Clayton (19871 ). 



3.3. A continuous low star formation rate 



As discussed by |Legrand et al. (1999|) , metals observed 
in IZw 18, and also in other starbursts galaxies, result 
from a previous star formation episode. Assuming that 
the present burst in IZw 18 is the first one, we evaluated 
the continuous star formation rate required to reproduce 
the observed oxygen abundance after 14 Gyr. In this sce- 
nario, a mild star formation process would have started a 
long time ago but the galaxy would presently undergo its 
first strong starburst event. We found that a SFR of only 
10 -4 g M Q yr" 1 (i.e., 10 -3 g M & Gyr^ 1 by unity of mass 
of gas in the galaxy), where g is the fraction of gas (in 
mass) available, can reproduce the observed oxygen abun- 
dance in IZw 18. Moreover, this model reproduces per- 



fectly the carbon abundance measured by Garnett et al 
[T997| ). The kinetic energy inject ion rate, evaluated us- 
ing the models of Cervino (1998| ), is for this scenario of 



910 ergs , i.e., probably insufficient to eject the met- 
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als out of the galaxy (Mac Low 
Tenorifr-Tagle 1998|). 



Fcrrara 1999: Bilich & 



Moreover, the kinetic energy is not deposited in one 
single region like in a burst, but as the continuous star for- 
mation is supposed to occur sporadically in location, the 
injected energy is diluted over the whole galaxy, reducing 



the efficiency of ejection of the metals (see also Strickland 



& Stevens 199E). For these reasons, we assumed that all 
the metals produced by the continuous star formation rate 
are retained in the galaxy. Of course, this continuous star 
formation regime represent an extreme case. We cannot 
rules out the existence of intermediate models in which 
the star formation history would be a succession of very 
small bursts without or with very low metal loss. However, 
these models will appear, in average, as a rather contin- 
uous star formation rate. The motivation for preferring a 
continuous star formation rate is principally the absence 
of observational evidences for gas rich galaxies with a SFR 
equal to zero, even among LSBG. 

Finally, in order to compare the colors predicted 
by the model, we added the present burst at 14 Gyr. 
The characteristics for the current burst, i.e., a SFR 
of 0.023 M yr _1 during 20 Myrs, were taken from 
Mas-Hesse fc Kunth (1999| ). 



were 

The evolution with time 



of the gas fraction, oxygen and carbon abundances is 
presented in Fig. pi whereas the evolution of the colors is 
shown in Fig. 0. 



Note that the observations of Thuan (1983) were done 
using a 8" circular aperture, which is smaller than IZw 18. 
Thus the total asymptoti c magnitudes may be smaller 
than the ones measured by Thuan (1983 ). Doublier (1998 ) 
has shown that this difference can be as large as 3 mag- 
nitudes, due to the presence of an old underlying stellar 
population. However, as IZw 18 is a very unevolved ob- 
ject, the old underlying population must be very faint. We 
have evaluated the difference between the observations of 



Thuan (1983 ) and the total magnitude expected from our 
model in the case of the existence of an old underlying 
population, due to continuous star formation, extending 
uniformly over the whole galaxy. We assumed that the old 
underlying population extends over 60" x 45" as discussed 
later. If m is the total magnitude emitted in a band, it can 
be written as 



m = -2.5Log(E b + E ci8 + E 



co8) 



(2) 



with Eb the flux emitted by the burst (localized in the cen- 
tral region included in the measurement of Thuan 1983| ), 
E c i8 and E co8 the fluxes emitted by the old underlying 
stellar population inside and outside the 8" aperture, re- 



spectively. The magnitude measured by Thuan (1983) are 
then: 



m 8 = -2.5Log(E b + E ci8 ) 



(3) 



E co8 can be evaluated using the flux of the old 
underlying stellar population predicted by the model. 
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Fig. 3. Time evolution of the gas fraction, oxygen and 
carbon abundances for the continuous SFR model with 
a Salpeter IMF. The dots represent the measured abun- 
dances. 
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Fig. 4. Colors evolution (U-B, B-V, V-K) for the contin- 
uous SFR model with a Salpeter IMF. The dashed lines 
delimit the zone of compatibility between the model and 
the observations. 



Under these assumptions, the magnitude measured by 
Thuan (1983 ) should be decreased by 0.20 in J and 0.24 
in K. Thus this does not change the main results and 
the model predicts colors consistent with the observations. 



The fraction of gas consumed remains very low; thus 
g is always close to 1. This means that the SFR is rather 
constant. The model is thus fully compatible with all the 
observations within the error bars. 



3. 4- First consequences 

If we suppose that the continuous star formation occurs 



sporadically over the whole galaxy like in LSBG (Van Zee 
| ct al. 1997c ), the observed homogeneity of abundances 



(within the NW region and also between NW and SE 
regions of IZw 18) is a natural consequence of this 
scenario. The uniformly distributed star formation and 
the long time evolution (14 Gyr) ensure the dispersal and 
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homogenizing of the metals over the whole galaxy. 

We evaluated the number of stars with a mass greater 
than 8 M Q , formed over 14 Gyr, to be around 12000. 
This corresponds to 120 massive stars (typically an open 
cluster) formed every 140 Myr. Taking their lifetime into 
account, we expect to see around 13 stars with mass 
greater than 8 M at a given epoch. We also evaluated 
the SN rate expected to be 7.5 10 -7 yr -1 , or relatively 
to the mass of gas, around 10~ 14 yr^.Mg 1 . This can be 
compared to the SN rate in our gala xy which amount 
10 -13 yr i M -i ( (Tammann et al. 1994|) . 



of the absorption measured by Kunth ct al. (1994[ ) is 
concentrated at density higher than 10 20 atoms cm^, we 
assumed that this level represent the density limit for 
the continuous star formation. This appears as a lower 
limit because metal produced by star formation in more 
inner regions can also have been dispersed and mixed at 
larger distances. This implies that local fluctuations in 
the density should be up to a factor of 150 and that the 
continuous SFR could occur over a surface of 60 x 45" 



( |Van Zee et al. 1998| ). 



3. 5. Threshold and efficiency of star formation 

There is a relationship between the star formation rate and 



The continuous star formation rate evaluated is very 
low. From the work of Wyse fc Silk (1989 ), we can compute 
the star formation efficiency of such a process. According 
to these authors, the star formation rate at a distance r 
of the center and a time t is: 



the gas surface density (Schmidt 1959). A simple power 
law describes this link at "high" density, but this relation- 



(6) 



ship breaks down under a critical threshold (Kcnnicutt 



1989, 1998). The existence of a gas density critical thresh- 



old under Bthich star formation wmilH hp inhihitpH haH 



where e is the star formation efficiency, f2(r) the local 
angular frequency and fijji (r, t) the surface d ensity in 
HI. Using the rotation curves computed by Van Zee 



been p "opospH hy [Quirk (1079) This threshold would hp 



associated with largp scalp gravitational insta.hilit.ips for 



iai 

it! 



et al. (199S) and Petrosian et al. (1997), we found that 
n = (8.76 10 7 yr) -1 . The HI mass measured by Van Zee 



formation nf massivp clouds Tn a thin isothermal dislr 



thc critical gas surface density threshold ( Toomrc 1964 ; 
Cowic 198l|) is flKcnnicutt 1989|) : 



3.36G 



with 

« = 1.41^(1+^/3 
R y VdR J 



(4) 



(5) 



where c is the dispersion velocity in the gas, k is the 
epicyclic frequency (derived from the rotation curve), and 
V the rotation velocity at the distance R fro m the center 
of the galaxy. However, Van Zee et al. (1997c ) have shown 
that despite a gas density lower than the threshold, 
LSBG are undergoing star formation. We suggest that 
in low density objects, the density can fluctuate locally 
and get above the threshold in some places. This could 
induce localized and faint star formation (Van Zee et al 



Zee et al. 



1997c; gkillman 1999). For IZw 18, observations of Van 



and Petrosian et al. (1997) reveal a 



solid body rotation (in the central part) with parameter 
^ = 70 kms^kpc -1 . Using a dispersion velocity in 
the gas of 12 kms -1 (Van Zee et al. 1998), the critical 



threshold is of the order of 1.5 10 22 atoms cm -2 . On the 
other hand, Van Zee et al. (1998| ) have shown that the 
abundances in the HI halo are comparable to the one in 
the HII region. This suggest that star formation may have 
also occurred quite far away from the central regions. 
This is reinforced by the observation of star form ation 
in regions with density lower than the threshold ( Van| 
Zee et al. 1997c). As most of the HI gas responsible 



et al. (199S) using a surface of 2.3x3 kpc gives a mean 
surface density of 3 M Q pc -2 . The star formation rate 
is then if) = 0.18 e M Q yr -1 , which compared with the 
results of our model (10 -4 Mq yr -1 ), leads to a very low 
star formation efficiency of e ~ 6 10~ 4 . 



4. Generalization of the continuous star 
formation hypothesis 

4-1. Underlying population and faint objects population 

If the starburst in IZw 18 is the first one in its history, 
we can expect that objects which have not yet undergone 
a burst (but only a continuous star formation rate) do 
exist. They should then look like IZw 18 just before 
the present burst. Our modeling predicts that after 
14 Gyr of continuous star formation, the magnitude of a 
IZw 18-like object is of the order of 20 in V and 17.5 in 
K. These magnitudes represent the brightness of the old 
underlying stellar population in IZw 18. Assuming that 
the continuous star formation occured in regions where 
the HI column density was greater than 10 20 cm -2 , i.e., 
60 x 45" ( Van Zee et al. 1998| ), the expected surface 
brightness would be of the order of 28 magarcsec 2 in V 
and 26 magarcsec 2 in K. These values are an upper limit 
(in magarcsec 2 ); if a fraction of metals is ejected out 
of the galaxy, the SFR needed to produce the observed 
abundances will be higher and the total luminosity 
and surface brightness will be increased. Moreover, as 



discussed in section 3.5, the density limit adopted for the 



continuous SFR is a lower limit and the region where the 
continuous SFR can occur may be smaller, resulting in 
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higher surface brightness. However, the extreme faintness 
of the old underlying population probably explains why 
no strong evidence for its existence has been found in 
IZw 18 ( [Thuan 1983| ; [Hunter fe Thronson 1995| ) until 
recently when reanalyzing HST archive images [Aloisi 
et al. ( |L99£ ) found stars older than 1 Gyr. In all the cases, 
these very low surface brightness levels will be reachable 
with 8m class telescopes like the VLT and Gemini for 
rather closeby objects. A search for a faint old underlying 
population, in the external part of dwarf galaxies like 
IZw 18, resulting from a continuous star formation 
process, is planned. If successful, this will support the 
existence of a class of very low surface brightness galaxies, 
which never underwent a burst, but evolved through a 
continuous weak star formation rate. 



4-2. LSBG as quiescent counterparts of starbursts 

We also studied the state of dwarf galaxies between bursts. 
IZw 18 appears as an extreme object which could present 
for the first time a strong star formation event. The chem- 
ical abundance levels in most of the starbursts galaxies 
suggest that they have undergone at least two or three 
previous bursts. Using the model described below, we in- 
vestigated the simple following star formation history for 
a dwarf galaxy with mass, size and distance comparable 
to IZw 18: 

— A continuous star formation rate since 14 Gyr. 

— Two bursts with a SFR of 0.065 M Q yr _1 and a dura- 
tion of 20 Myr, respectively at 8 and 11 Gyr. 

The absolute magnitudes predicted are around -10 and 
the surface brightness about 25 magarcsec -2 in B. These 
values correspond to objects at the extreme en d of the 
luminosity function of the galaxies as observed by Lo et al 



(1993| ) 



We have also compared the continuous SFR required 
for IZw 18 (10~ 4 Mq yr^ 1 ) with those measured by Van 
Zee ct al. (19974 p]R) in low surface brightness galaxies 



As these objects have different sizes and masses, we nor- 
malized the SFR to the total HI mass. For IZw 18, the HI 



mass lies between 2.6 10 7 M ( [Lequeux fc Viallcfond 1980| ; 
Van Zee ct al. 1998|, for the main body) and 6.9 10 7 M© 
( Van Zee ct al. 1998| , for the whole galaxy, including the 
diffuse low column density component). The comparison 
is shown in table [l]. It appears that the continuous SFR 
as predicted by our scenario is comparable, relative to the 
HI mass, to the lowest SFR observed in quiescent and low 
surface brightness galaxies (like for example in UGC8024 
or UGC9218). 

We then conclude that LSBG and quiescent dwarfs 
are lik ely to be the quiescent counterparts of starburst 
galaxies, 



Table 1. Comparison between the SFR measured in quies- 
cent dwarfs and LSBG and the continuous SFR predicted 
for IZw 18. 



Name 


T ^ a 

ypc 


Mb 


iVi HI 
10 7 M R 


SFR 
M /yr 




SFR 

M HI 


rcf b 


UGCA20 


Irr/LSBG 


-14.9 


25 


0.0083 


3 




11 


V97c 


UGC2684 


dW/LSBG 


-13.7 


15 


0.0015 


1 


OOe 


11 


V97c 


UGC2984 


7/LSBG 


-18.4 


500 


0.35 


7.00c 


11 


V97c 


UGC3174 


Irr/LSBG 


-15.7 


61 


0.0086 


1 


41e 


11 


V97c 


UGC5716 


dSp/LSBG 


-16.3 


140 


0.023 


1 


64e 


11 


V97c 


UGC7178 


dl/LSBG 


-16.6 


140 


0.022 


1 


57e 


11 


V97c 


UGC11820 


dSp/LSBG 


-17.7 


440 


0.074 


1 


68e 


11 


V97c 


UGC191 


dSp/ qui 


-18.2 


320 


0.15 


4 


69e 


11 


V97c 


UGC300 


dwf/LSBG 


-16.5 


93.3 


0.03 


3 


26e 


11 


V97b 


UGC521 


Irr/LSBG 


-15.8 


95.5 


0.023 


2 


41e 


11 


V97b 


UGC634 


dSp/ qui 


-17.7 


470 


0.029 


6 


17e 


12 


V97c 


UGC891 


dWf/qui 


-16.3 


100 


0.014 


1 


40e 


11 


V97c 


UGC1175 


?/qui 


7 


102 


0.0037 


3 


63e 


12 


V97b 


UGC2162 


dl/qui 


? 


71 


0.0059 


8 


31e 


12 


V97b 


UGC2535 


dl/qui 


? 


275 


<0.003 


1 


09e 


12 


V97b 


UGC3050 


S/qui 


? 


479 


0.13 


2 


71e 


11 


V97b 


UGC3672 


Irr/LSBG 


-16.6 


129 


0.021 


1 


63e 


11 


V97b 


UGC4660 


dSp/LSBG 


-17.8 


363 


0.025 


6 


89e 


12 


V97b 


UGC4762 


dl/qui 


? 


100 


0.0099 


9 


90e 


12 


V97b 


UGC5764 


dl/qui 


-13.9 


34 


0.0098 


2 


88e 


11 


V97c 


UGC5829 


dl/qui 


? 


182 


0.080 


1. iii< 


11 


V97b 


UGC7300 


dl/qui 


? 


129 


0.016 


1 


24e 


11 


V97b 


UGC8024 


dl/qui 


? 


49 


0.0019 


3 


88e 


12 


V97b 


UGCA357 


dl/LSBG 


-16.3 


190.5 


0.034 


1 


79e 


11 


V97b 


UGC9128 


dl/qui 


? 


3.55 


0.00017 


4 


79e 


12 


V97b 


Haro43 


S /LSBG 


-17.1 


2 20 


0.092 


4 


02e 


11 


V97b 


UGC9762 


dSp/LSBG 


-18.0 


479 


0.05 


1 


04e 


11 


V97b 


UGC10281 


dwf/LSBG 


-17.1 


67.6 


0.012 


1 


78e 


11 


V97b 


IZw 18 c 
IZw 18 c 


dl 
dl 


-13.92 
-13.92 


6.9 
2.6 


0.0001 
0.0001 


1 

3 


45c 
85c 


12 
12 





ological type / LSBG or quies cent dwarf fcrni^ 



a Morpl 

6 V97b: Van Zee et al. (1997b| ); V97c: |Van Zee ct al. (1997c ) 

c The two values given depends if the total or only the main component 

HI mass is adopted. 



4-3. Generalization to all galaxies 

If a continuous star formation rate exists in IZw 18, it must 
exist in other dwarf galaxies, and may be, in all galaxies. 
There are some hints for such a hypothesis. 

For example, the extreme outerparts of spirals, where 
no strong star formation event occurred and where the 
metals formed in the more active inner zones have not dif- 
fused, must have low abundances and low surface bright- 
ness. For example, let us assume that the bulk of "strong" 
star formation occurs in a spiral galaxy at distances less 
than the optical radius (~10 kpc). Roy & Kunth (1995) 
have shown that metals can be dispersed at scales up to 
10 kpc in about 1 Gyr, so extending their results we can 
expect that if recent (less than 1 Gyr ago) star formation 
occurred in external regions located at 10 kpc from the 
center (one optical radius), the newly formed metals could 
affect abundances at distances up to 20 kpc form the cen- 
ter in few Gyrs. If this star formation is relatively recent, 
we expect that the most external region of the disk (more 
than 2-3 optical radii) will not be affected by "strong" star 
formation events their mctallicity will be solely the result 
of the "underlying" low continuous star formation rate. 
Extrapolations of metallicity gradients in spiral galaxies 
lead to abundances comparable to that of IZw 18 at ra- 



dial distances of about three optical radii ( Ferguson et al 
1998; Henry fc Worthey 1999| ). This corresponds to the 



size of the halos or disks suscept ible to give rise to meta l 
lie absorption in quasar spectra (Bergeron fc Boisse 1991) 
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. Prochaska & Wolfe (99) 




A Lu et al (96) 











2 3 

Redshifl z 



Fig. 5. Comparison of the predicted and observed evolu- 
tion with redshift of the abundance [Fe/H]. The points 
represent the data from Lu ct al. (1996) and prochaska fc 
Wolfe 1999) and the solid line the model prediction for a 
constant star formation rate. 



We showed that a continuous low star formation rate 
results in a steady increase of the metallicity of the inter- 
stellar gas. We have compared the evolution of the iron 
abundance predicted from our model with the measure- 
ments in DLA systems in Fig. 0. The abundances pre- 
dicted by the model mimic the lower envelope of these 
measurements. If we assume that these absorption sys- 



tems are associated with galaxy halos ( Lanzctta ct al 
1995; fcipp et al. 1997 ), this indicates that such a process 
can account for a minimal enrichment of the ISM. One 
measurement appears lower than the model prediction. 
However, Fe atoms arc likely to condensate into grains 
( Lu et al. 1998| ; Pettini ct al. 1997), so the iron abundance 
measurements are only lower limits of the real iron abun- 



dances. Moreover Bergeron & Boisse (1991) have shown 



that absorptions should occur at distances of up to 4 
Holmberg radius. These regions are likely to present very 
low density, and may be too under-critical to allow star 
formation, even in the low SFR regime described here. 
Moreover, the partial ioniza tion of these regions by the 
diffuse ioni z ing background ( Van Gorkom 1991 ; Corbclli 



et al. {98|; |Maloney 1990| ; |Corbelli fc Salpeter 1993a| ]b|) 



will also contribute to prevent star formation. Their metal- 
licity could thus be due only to metals which have diffused 
from the inner regions. If true, we can expect, in these re- 
gions, abundances lower than what is predicted from the 
continuous SFR. 



5. Conclusion 

We have investigated different star formation histories for 
IZw 18 using a spectrophotometric model coupled with 
a chemical evolution model of galaxies. We have shown 
that if the observed metallicity results only from burst 
events with galactic winds, no more than 50-70 % of the 
newly synthesized metals may have been ejected out of the 
galaxy. This is because a larger metal loss rate will require 
to form more stars to reach the measured abundances, 
hence resulting in redder colors than what is observed, 
due to an overproduction of old underlying low mass stars. 
Following the suggestion of Legrand et al. (1999), we in- 
vestigated the hypothesis of a low, but continuous, SFR 
which should account alone for the observed metallicity 
in IZw 18. We have shown that the metals in IZw 18 are 
likely to result from a mild continuous star formation rate 
which took place indenpently from bursts. This star for- 
mation would be due to local fluctuations in the density 
which exceeds sporadically the threshold for star forma- 
tion. Using a spectrophotometric model and a chemical 
evolution model of galaxies, we demonstrated that a con- 
tinuous star formation rate as low as 10 -4 M^lyr occur- 
ring for 14 Gyrs can reproduce all the main parameters 
of IZw 18. The generalization of this model to all galax- 
ies accounts for many observed facts, such as the presence 
of star formation in quiescent dwarfs and LSBG, the in- 
crease with time of the metallicity of the most underabun- 
dant DLA systems, the extrapolation of metallicity in the 
outerparts of spiral galaxies, the lack of galaxies with a 
metallicity lower than IZw 18, the apparent absence of 
HI clouds without optical counterparts, and the homo- 
geneity of abundances in dwarfs galaxies. Moreover, we 
predict for IZw 18 and other extremely unevolved galax- 
ies, the presence of an old underlying stellar population 
(resulting from this continuous star formation process) at 
a surface brightness level of at least 28 mag arcsec 2 in V 
and 26 mag arcsec 2 in K. Finally, we have shown that the 
parameters for the low continuous star formation rate are 
comparable to what is observed in LSBG and quiescent 
dwarfs, suggesting that these objects could be the quies- 
cent counterparts of starburst galaxies. 
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